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Abstract 

QCD corrections to the penguin induced charmless inclusive B decays are performed 
to an improved leading-log approximation. For process b — * sqq (q=uds) and bq — > sq, 
four previous missed 4-quark operators are included to make a complete leading log QCD 
corrected result. Furthermore, part of the next-leading log effect was also given for the above 
processes and process b — ► sg. In comparison to previous QCD corrected calculations by 
Grigjanis et al., the inclusive decay rates for process b — > sqq (q=ud), b — * sss and bq — > sq 
are suppressed 10%, 14% and 18% respectively. The branching ratio of process b — > sg 
is enhanced by 2-3 orders, the result for this process is in agreement with the result of 
Ciuchini et al. 
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1 Introduction 



Recently the CLEO collaboration has observed the charmless decays of B, B — > 7r + 7r~, B° — > 
K + ir~ and B° — > The upper limits on the branching fractions are f?^ < 2.9 x 10 -5 , 

Bk-k < 2.6 x 10~ 5 , Bkk < 0.7 x 10~ 5 . The sum of B nn and Bx n exceeds zero with a significance 
more than four standard deviations. One can expect, individuals of these exclusive channels and 
also the inclusive processes will be found at CESR or future B factories. 

The B meson system is expected to give rise to a very rich phenomenology, providing a wealth 
of information that will allow us to constrain the currently very successful standard model (SM), 
as well as many extensions that go beyond it. It also provides many channels for the search of 
CP violation. 

The charmless b quark decay has already been calculated in several papers @, ||, £|, ||, ||. They 
are induced by both b — > u transition and penguin diagrams. The b — > u transition diagrams 
are easy to calculate, but penguin diagrams are more complicated. W.S. Hou showed the results 
of penguin diagrams without QCD corrections in ref.0. Afterwards, results including some 
short distance QCD corrections were given by R. Grigjanis et al.Q. For different processes, the 
authors claimed that the branching ratios are either QCD-enhanced or suppressed by a notable 
factor. In other words, the strong interaction plays an important role in these decays. But this 
is not the end of the story, there are still some unsolved questions. The anomalous dimensions 
used in ref. H were changed later by the authors themselves in ref.0's erratum. This change 
caused large effects in process b — > sg which can be seen in the present paper. See also a recent 
work by Ciuchini et al. JS|. Furthermore, only two four-quark operators were considered there, 
while the other four were neglected. Unfortunately, the missing four operators are responsible for 
contributing to the matrix elements of processes b —>■ sqq (q=uds) and bq — > sq. Although their 
coefficients of matching at Mw can be transferred to the operator O7 of ref. HQ, when running, 
these six four-quark operators will be mixed together, so the missing four will reappear with 
3 Note this operator O7 is not the one that we defined in this paper 
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nonzero coefficients. This is similar to the b — > 57 case.Q 

Since charmless decays of B meson caused a great interest in studying CP violation, and 
the present experimental limit falls in with the theoretical predictions of the standard model, a 
more accurate calculation is needed to reduce theoretical uncertainties. The aim of the present 
paper is to present an improved leading log QCD corrected result. 



2 Matching at M } 



w 



The basic effective field theory idea is by now quite well established |T0|, pjj . In our particular 
case of Minimal Standard Model, we integrate out the top quark and the weak W bosons at 
jj, = Mw scale, generating an effective five-quark theory. By using the renormalization group 
equation, we run the effective field theory down to b-quark scale to calculate QCD corrections 
to charmless b decay. 

After applying the full QCD equations of motionfQ, a complete set of operators relevant 
for charmless B decays can be chosen to be: 

0\ = (c L ^b La )(s La -f^c L p), 2 = (c La Yb La )(s L ^^c L p), 

q q 

q q 
7 = {gz/lQK 2 )m b -s L a^T a b R Gl v . 

The covariant derivative is defined as 

with g 3 denoting the QCD coupling constant. 

Then we can write down our effective Hamiltonian as 

H ef f = 2V2G F V tb V* ]T CiMOiQi). (2) 



4 See discussions in ref . |fj, |j| . 



One can find the coefficients of operators at fi = Mw scale by integrating out the weak gauge 
bosons and would-be Goldstone bosons at this scale. Since in standard model, there is no 
tree level flavor changing neutral current, all charmless b — > s processes are through loop dia- 
grams. They are QCD induced processes. Even if no radiative QCD corrections are performed, 
they are already order of 0(a s )@. The decay amplitudes of four-quark processes are propor- 
tional to [f\a s \og{Mw / TTib) 2 + /2«s]- The leading logarithmic QCD corrections to them are of 
0([a s log(Mw /iTLb) 2 ] n ) with n > 2. They are summed by the one-loop renormalization group 
equation in the leading log approximation. So this approximation will include the leading log 
terms of all loops but miss the f2Q. s terms of one-loop diagrams. Although this is usually con- 
sidered as part of the next-to- leading log effect, it is actually part of the one- loop result. We 
will use the complete one-loop result (including f'2<y s ) for four-quark operators in matching at 
M\y, so that we can get a result from both the complete one- loop contribution and summation 
of leading log contribution. This is an improvement in the perturbative calculation at least for 
the QCD-induced four-quark processes. The matching diagrams for four-quark operators are 
displayed in Fig.[TJ (There should also be three additional symmetric four quark diagrams at each 
side of the first equation in Fig.l.). Matching diagrams for operator Oj are very similar to b —>■ 57 
calculations 0]. Neglecting small terms proportional to m 2 or rn? u in these matching conditions, 
and using V c bV* s = —V t bV t * s , one finds the following coefficients of the operators: || [TJ 



C 1{ M W ) = Ua j Mw \ C 2 (M W ) = x^WM.) 



in ' ~^'"™> 24tt ' 

C 3 (M W ) = C 5 (M W ) = ^^-F X {5) 

C,{M W ) = C 6 {M W ) = -^M^l Fl (6) 



where 



™=f + v^ 2 + (i -V iog * (4) 



with 5 = M^/m 2 . 
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Table 1: Numerical results for coefficients of operators Ci(m&) with a s {mz) = 0.117. 
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3 Renormalization group running from Mw to 



We then use the renormalization group equation satisfied by the coefficient functions Ci(fi), to 
continue running the coefficients of operators from /i = My/ to \i — m^. 



(5) 



Where %jS are anomalous dimensions of operators. These anomalous dimensions have been 
calculated by many authors for the process b — > 37. Only recently it is completely solved by 
Ciuchini et al.|15|], and their calculation is confirmed by Cella et al.[|T6|j. 
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In the leading order of $3, the solution to eqn.([5[) in matrix notation is given by 



C(/i 2 



9 M 7 r 

exp / d 9-TTT 



CO* 



(6) 



(7) 
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After insertion of anomalous dimension matrix(^), we get the coefficients of operators at \i — m b 
scale. Here we made use of Mw = 80.22GeV, m b = 4.9GeV, a s {mz) = 0.117[|T7|]. Table [L] shows 
the numerical values of coefficients of operators Oi, with different input of top quark mass. In this 
table, coefficients of four-quark operators 3 ,4, 5 , 6 have very little dependence on the top quark 
mass, while C 7 (m b ) is getting bigger when top mass increases. C^fm?,) = 1.077. Its coefficient 
does not vary with top mass at all, because operator O2 is generated from W exchange diagram 
and does not mix with top quark loop induced operators, 



4 The charmless b decay rate 



There are four types of nonleptonic charmless decays considered here. Their inclusive decay 
widths are given by the sum of operators Oi in our effective field theory. 



(1) Operator O5 and O7 in our operator basis contribute to process b — > The decay 

width is, 

8a, 



7T 



c 7 e// K) 



(8) 



where 



1 G 2 F m b 2 



'ts v tb\ , 



192vr 3 

C e 7 ff (m b ) = C 7 (m b ) + C 5 (m b ). 

(2) For process b — ► sqq(q denotes u or d quark), the Feynman diagrams are displayed in 
Fig.^. Operators 02,3,4,5,6,7 contribute to this process in the effective field theory. The decay rate 
is. 



*sqq 



3]T \cf f (m b )\ 2 + 2 {C e 3 ff (m b )C e / f (m b ) + C e 5 ff (m b )C e 6 ff (m b )} 



i=3 



r , 



(9) 



where 



Cf f (m b ) = d 



a* , m 



m b ) 



log— § C 2 (m b ), i = 3,5, 



36tt ml 
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Cf f {m h ) = Ci{m b ) + T^log^f C 2 (m b ), % = 4,6 



a 



12tt 



with m c = l.hGeV . Here we included all terms proportional to a s , to make a complete one-loop 
0(a s ) result. 

(3) For process b — ► sss, there are additional diagrams concerning momentum exchange of 
two s quarks. The decay width is, 



b—>sss 



Cl ff (m b ) + Ct ff {m b ) 2 + 3 C e 5 ff (m b ) 2 + 3 C e 6 ff (m b ) 



(10) 



+2C e 5 ff (m b )C e e ff (m b ) + {C e / f (m b ) + C e / f (m b ) + C e 6 ff (m b )}C 7 (m b ) 



rv 



(4) The Feynman diagrams contributing to the nonspectator process bq — > sq are also seen 
in Fig.||. The decay width is, 

2 



r 



bq^sq 



327^ 



c 5 e// K) + 3C 6 e// K) 2 (^] r . 



Where the B decay constant is fs = 200MeV. 

In order to find the branching ratios of these kinds of charmless B decays, the semileptonic 
decay of B is used. 

rfl 



BR(B — > X s no charm) 



Sn ° Charm) -BR(B^X c eV) 



T{b —> ceu) 



where 



r(6 -> cei/) ~ g(m c /m b ) M - ^^/K/m 6 )j r , 



(12) 



(13) 



with g(m c /m b ) ~ 0.447 and f(m c /m b ) ~ 2.4 correspond to the phase space factor and the 
one-loop QCD correction to the semileptonic decay, respectively [113]. Here we use experimental 



result J Br(5 -> X c eT7) = 11% |7 



The branching ratios for different processes are given in Fig. In comparison to the previous 
QCD-corrected results of ref.||, the decay rate of process b — * sg is strongly enhanced by QCD 
corrections rather than severely suppressed. The differences are of 2-3 orders. The reason for such 
a large difference, is obvious, for the authors of ref.p| used the wrong anomalous dimensions of 



ref.JTJ, which has already been corrected in its erratum. If right values are taken, the differences 
are not so large. Our result for b — > sg is consistent with the recent work by Ciuchini et al.|| 
in leading log approximation. The part of next-to-leading log contributions we considered in 
this paper is not essential to this process. For process b — > sqq, b — > sss and bq — > sq, after we 
include all the dimension-6 four-quark operators, the complete leading log QCD corrections make 
the branching ratios suppressed by 15%, 19%, and 25% respectively (comparing to the result of 
ref-H), when m tov = 175GeV. Furthermore, we also include 0(a s ) terms in the matching 
condition, which is the first term of next-to-leading log effects. The branching ratios for these 
processes are slightly enhanced. Comparing with ref.||, the total effects are 10%, 14%, and 18% 
suppression of the branching ratios, respectively, 

The branching ratios obtained by W.S. Hou|2[ without QCD running from Mw to are 
also given at Fig. |J Comparison with this result shows that our decay rate of process b — > sg 
is strongly enhanced, e.g. a factor of 3.4 at m top = 175GeV, this is almost the same as b — ► 57 
decay fTB"|, |l~5f . While for process b — > sqq, b —>■ sss and bq — > sq, the branching ratios are only 
increased by 49%, 34% and 17%, respectively. 

5 Conclusion 

In conclusion, we have given the full leading log QCD corrections to penguin induced charmless 
b —>■ s processes together with some next-leading log corrections. The branching ratios of pro- 
cesses b — > sqq (q=uds), and bq — > sq are found to change slowly with the top quark mass. For 
process b — > sg, the decay rate is getting larger when top mass increases. The result also shows 
that for b —>■ sg process(gluon on shell), QCD corrections are as important as in the b — > s'j case. 
The QCD corrected results for four-quark processes are all enhanced, the amount of the increase 
vary notably for different processes. 
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Figure Captions 

Fig. 1: Matching conditions at /i — M w for Green functions in the full standard model(left 
hand side) and effective field theory below W scale (right hand side) with the heavy 
dots denoting high dimension operators. 

Fig. 2: Feynman diagrams contributing to b — > sqq, bq — > sq in effective field theory. 

Fig. 3: QCD corrected branching ratios as function of top quark mass. From top to bottom 
are lines for process b — > sg, b — > sqq, b — > sss and bq sq. 

Fig. 4: Branching ratios without QCD radiative corrections, obtained by W.S. Hou as function 
of top quark mass. From top to bottom are lines for process b — > sqq, b — > sss, bq — ► sq 
and b — > sg. 
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